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ABSTRACT 

We present observations of the Sunyaev-Zel'dovich effect (SZE) in the Bullet cluster (IE 
0657-56) using the APEX-SZ instrument at 150 GHz with a resolution of 1'. The main results 
arc maps of the SZE in this massive, merging galaxy cluster. The cluster is detected with 23 a 
significance within the central 1' radius of the source position. The SZE map has a broadly 
similar morphology to that in existing X-ray maps of this system, and we find no evidence for 
significant contamination of the SZE emission by radio or IR sources. In order to make simple 
quantitative comparisons with cluster gas models derived from X-ray observations, we fit our 
data to an isothermal elliptical /3 model, despite the inadequacy of such a model for this complex 
merging system. With an X-ray derived prior on the power-law index, /3 = 1.041;q }q, we find a 
core radius Tc = 142 ± 18", an axial ratio of 0.889 ± 0.072, and a central temperature decrement 
of —771 ± 71/iKcMB, including a ±5.5% flux calibration uncertainty. Combining the APEX- 
SZ map with a map of projected electron surface density from Chandra X-ray observations, 
we determine the mass-weighted temperature of the cluster gas to be Tmg = 10.8 ± 0.9 keV, 
significantly lower than some previously reported X-ray spectroscopic temperatures. Under the 
assumption of an isothermal cluster gas distribution in hydrostatic equilibrium, we compute the 
gas mass fraction for prolate and oblate spheroidal geometries and find it to be consistent with 
previous results from X-ray and weak lensing observations. This work is the first result from the 
APEX-SZ experiment, and represents the first reported scientific result from observations with 
a large array of multiplexed superconducting transition-edge sensor bolometers. 

Subject headings: cosmic microwave background — cosmologyiobservations — galaxies: clusters: indi- 
vidual (IE 0657-56) 

1. Introduction 
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In particular, active mergers of sub-clusters pro- 
vide a window to the processes by which mas- 
sive clusters are assembled. In these systems, the 
galaxies and associated dark matter are essentially 
coUisionless. In contrast, the ionized intrac luster 
gas, typically at temperatures of T ~ 10^ K, is 
strongly interacting and experiences complex dy- 
namics. In extreme cases, the normally associated 
dark matter and intracluster gas can be signifi- 
cantly separated. 

The Bullet cluster (IE 0657-56) at z = 0.296, is 
a massive cluster consisting of two sub-clusters in 
the process of merging. The smaller sub-cluster or 
"bullet" has passed through the larger main clus- 
ter. X-ray observations infer a bow shock velocity 
of ~4700 km/s (Markevitch 2006), while simula- 
tions of the collision yield a substantially lower 
speed for the sub-cluster (Springel & Farrar 2007). 
This collision is perpendicular to the line of sight, 
providing an ideal system for studying interacting 
sub-clusters (Clowe et al. 2006). 

The mass surface density of the Bullet cluster 
has been measured using weak and strong gravi- 
tational lensing of light from background galaxies. 
There are significant angular offsets between the 
peaks of the X-ray surface brightness, which trace 
the baryonic gas through thermal bremsstrahlung 
emission, and the peaks of the lensing surface den- 
sity, which are associated with the majority of 
the mass. The combined weak and strong lens- 
ing analyses of Bradac et al. (2006) show that the 
main cluster and sub-cluster are separated from 
their associated X-ray peaks at 10 a and 6 a sig- 
nificance respectively. This result has been recog- 
nized as providing direct evidence for the presence 
of coUisionless dark matter in this system (Clowe 
et al. 2006). 

The Sunyaev-Zel'dovich effect (SZE) provides 
an independent probe of the intracluster gas. In 
the SZE, a small fraction (^ 1%) of cosmic mi- 
crowave background (CMB) photons undergo in- 
verse Compton scattering from intracluster elec- 
trons (Sunyaev & Zel'dovich 1970; Birkinshaw 
1999). This process distorts the Planck black- 
body spectrum of the CMB and produces a signal 
proportional to the gas pressure integrated along 
the line of sight. At 150 GHz, the SZE produces 
a temperature decrement with respect to the un- 
perturbed CMB intensity. Early detections of the 
SZE in the Bullet cluster include work by An- 



dreani et al. (1999) and Gomez et al. (2004). 

Unlike the X-ray surface brightness, the peak 
SZE surface brightness for a given cluster is in- 
dependent of redshift. Therefore, the SZE has 
the potential to be an effective probe of intra- 
cluster gas out to the redshifts at which clusters 
are assembled. SZE measurements of galaxy clus- 
ters provide complementary constraints on cluster 
properties typically derived from X-ray measure- 
ments such as central electron density, core radius 
of the intracluster gas, cluster gas mass, and frac- 
tion of the total cluster mass in gas. Since the SZE 
and X-ray signals are proportional to the line-of- 
sight integral of the electron density and electron 
density squared respectively, SZE results will be 
less sensitive to clumping of the intracluster gas. 
For all comparisons between SZ and X-ray data, 
we assume a ACDM cosmology, with h = 0.7, 
f^,„ = 0.27, and f^A = 0.73. 

In this paper, we present a 1' resolution SZE 
image of the Bullet cluster at 150 GHz made with 
the APEX-SZ instrument. It is the first reported 
scientific result from observations with a large 
array of multiplexed superconducting transition- 
edge sensor bolometers. In §|2] we discuss the 
instrument and observations. Calibration is dis- 
cussed in §[3] In §|4] we describe the data reduction 
procedure, and in §[5) we present the results of fits 
to the SZE surface brightness with cluster mod- 
els, including mass- weighted electron temperature 
and gas mass fraction calculations. We summarize 
the conclusions and discuss future work in §[6j 

2. Observations 

APEX-SZ is a receiver designed specifically 
for SZE galaxy cluster surveys (Schwan et al. 
2003; Dobbs et al. 2006; Schwan et al. 2009, 
in preparation). It is mounted on the 12-meter 
diameter APEX telescope, located on the Ata- 
cama plateau in northern Chile (Giisten et al. 
2006). The observing site was chosen for its ex- 
tremely dry and stable atmospheric conditions. 
The mean atmospheric transmittance is frequently 
better than 95% in the APEX-SZ frequency band 
at 150 GHz (Peterson et al. 2003; Chamberlain 
& Bally 1995). The telescope is capable of round- 
the-clock observations. 

Three reimaging mirrors in the Cassegrain 
cabin couple the APEX telescope to the focal 
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plane of APEX-SZ. We achieve the diffraction 
Umited performance of the telescope across the 
entire 0.4° field of view with a mean measured 
beam full width half maximum (FWHM) of 58", 
and a measured beam solid angle of 1.5 arcmin^, 
including measured sidelobes at the —14 dB level. 

The APEX-SZ receiver houses a cryogenic focal 
plane, operating at 0.3 K. The focal plane contains 
330 horn-fed absorber-coupled superconducting 
transition-edge sensor bolometers (Richards 1994; 
Lee et al. 1996), with 55 detectors on each of six 
sub-array wafers. Of the 330 detectors, 280 are 
read out with the current frequency-domain multi- 
plexed readout hardware. We measure the median 
individual pixel Noise Equivalent Power (NEP) to 
be 10^^^ W/\/Hz and the median Noise Equiv- 
alent Temperature (NET) to be 860 /xKcmb\/s- 
The measured optical bandwidth of the receiver 
is 40% narrower than the design goal of 38 GHz, 
resulting in lower sensitivity than anticipated. 

The large field of view of the APEX-SZ instru- 
ment is designed for surveying large areas of sky. 
In order to efficiently observe a single target, we 
use the circular scan pattern illustrated in Fig- 
ure [l] The circle center is fixed in AZ /EL coor- 
dinates for twenty circular sub-scans, with a total 
duration of 100 seconds. This choice has a num- 
ber of important advantages. The sky signal is 
modulated so that it appears in the timestream 
at frequencies higher than atmospheric drifts and 
readout 1// noise. In addition, the circle scan has 
a moderate continuous acceleration; the lack of 
high acceleration turn-arounds makes it possible 
to achieve a high observing efficiency. Approxi- 
mately 20% of the total observing time is spent 
moving the telescope to a new center position be- 
fore the start of the next scan. Every bolometer 
maps a 12' x 25' sub-field, with a combined map 
field of 36' x 48' every 100 seconds. 

Observations of the Bullet cluster were con- 
ducted over a period of seven days in August 2007, 
when the cluster was visible between the hours of 
03:00 and 15:00 local time. The weather over this 
period was typical for the site, with precipitable 
water vapor varying between 0.25 and 1.5 mm, and 
a median atmospheric transmittance of 97%. For 
the analysis in this paper, 235 scans are used, each 
scan consisting of twenty 5-s circular sub-scans, for 
a total of 6.4 hours of on-source data. 
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Fig. 1. — The 100-second circular drift scan pat- 
tern. The solid line shows the track of the center of 
the array. One circle has a period of five seconds. 
The dashed line shows the instantaneous field of 
view of the bolometer array. The -I- marker indi- 
cates the source position with respect to the scan 
pattern. The center of the circles is constant in az- 
imuth and elevation as the source drifts across the 
field. The small disk in the lower left indicates the 
58" mean FWHM beam for a single bolometer. 

3. Calibration 

The response of the receiver to astronomical 
sources is measured with daily raster scans of 
Mars over every bolometer in the array. For each 
bolometer, the observations provide a primary flux 
calibration and a high signal-to-noise beam profile 
from which we determine beam parameters such 
as size, ellipticity, and position with respect to the 
array-center pointing. Additional observations of 
RCW57 and RCW38 are used to monitor gain sta- 
bility, and frequent observations of bright quasars 
near the cluster source are used to monitor point- 
ing stability. 

The WMAP satellite has been used to calibrate 
the brightness temperature of Mars at 93 GHz in 
five measurement periods spanning several years 
(Hill et al. 2009). The WMAP Mars tempera- 
tures are tied to the CMB dipole moment and are 
accurate to better than 1.0%. The brightness tem- 
perature of Mars changes significantly (~ 15%) as 
a function of its orbit and orientation. We use 
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a version of the Rudy Model (Rudy et al. 1987; 
Muhleman & Berge 1991), that has been updated 
and maintained by Bryan Butler to transfer the 
WMAP Mars temperature results to the APEX- 
SZ frequency band and specific times of our Mars 
observations. 

After applying a constant scaling factor, we find 
the Rudy model predictions for the Mars bright- 
ness temperature to be in excellent agreement with 
the WMAP measurements. We find that the Rudy 
model brightness temperatures at 93 GHz arc sys- 
tematically a factor of 1.052 ± 0.010 higher than 
those measured by WMAP in the five published 
observation periods. In contrast, repeating the 
same exercise with the thermal model developed 
by Wright (1976, 2007), as implemented in the on- 
Hne JCMT-FLUXES programQresults in a scaling 
factor of 1.085 ±0.043. This is consistent with the 
10% rescaling of this model called for in Hill ct al. 
(2009), but the scaling factor exhibits significantly 
larger rms scatter than that of the Rudy Model. 

We therefore use the WMAP 93 GHz cahbrated 
Rudy Model to compute the Mars brightness tem- 
peratures at 150 GHz for the specific times of our 
Mars observations by reducing the Rudy model 
150 GHz temperatures by a factor of 1.052. The 
Rudy model 93 GHz to 150 GHz frequency scaling 
factor is 1.016 ± 0.009 at the times of our Mars 
observations, and we adopt the rms scatter in this 
frequency scaling factor as an estimate of its uncer- 
tainty. Combining the uncertainties in the WMAP 
Mars calibration, the WMAP to Rudy model scal- 
ing factor at 93 GHz, and the Rudy model fre- 
quency scaling factor, we estimate the uncertainty 
in Mars temperature to be ±1.7%. 

The measured signals from the calibrators are 
corrected for atmospheric opacity, which is mea- 
sured with a sky dip observation at the begin- 
ning and end of each day's observations. Mea- 
sured zenith transmittance over the observing pe- 
riod ranged between 0.92 and 0.98, with a median 
of 0.97. Based on the observed temporal variabil- 
ity of the opacity, drifts in atmospheric opacity 
between the sky dip and observation contribute 
< 0.4% to the overall calibration uncertainty. Af- 
ter correcting for the atmospheric opacity, we find 
that the Mars temperature measured by APEX-SZ 
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varies from the model prediction by up to ~ 3% 
over the course of the observation period. This 
gain variation is included as a source of error in 
the final calibration uncertainty. The APEX-SZ 
observing band center is measured with a Fourier 
transform spectrometer to be 152 ± 2 GHz. The 
uncertainty in the band center results in a ±1.4% 
uncertainty in extrapolation of the Mars based cal- 
ibration to CMB temperature. 

The beam shape, including near sidelobes, is 
characterized by creating a beam map from Mars 
observations, combining the same bolometer chan- 
nels that are used to make the science maps. We 
adjust the calibration and measured beam size for 
the small 1%) correction due to the 8" angular 
size of Mars. We estimate a fractional uncertainty 
in the beam solid angle of ±4%. 

The APEX-SZ detectors operate in a state of 
strong negative electrothermal feedback which re- 
sults in a linear response to changes in the in- 
put optical power. We have measured the re- 
sponse of the detectors during sky dips between 
90° and 30° elevation (antenna temperature dif- 
ference ^ 13 K), and find no significant devia- 
tion from the expected linear response to loading. 
We therefore conclude that detector non-linearity 
makes a negligible contribution to the calibration 
uncertainty. 

Slowly changing errors in telescope pointing re- 
sult in both a pointing uncertainty and a flux cal- 
ibration uncertainty due a broadening of the ef- 
fective beam pattern. To measure pointing er- 
rors during our observations, we observe a bright 
quasar within a few degrees of the Bullet cluster 
every 1-2 hours, and apply a pointing correction 
as needed. The typical RMS pointing variations 
of the APEX telescope between quasar observa- 
tions is ~ 4". This pointing uncertainty results 
in a slightly larger effective beam for the coadded 
maps than is measured with the individual cali- 
brator maps. The correction to the flux calibra- 
tion of the coadded maps due to pointing uncer- 
tainty is negligible, particularly for the observation 
of extended objects such as the Bullet cluster. We 
estimate the pointing uncertainty in the coadded 
maps to be ±4" in both the RA and DEC direc- 
tions. 

The uncertainty in the CMB temperature cali- 
bration of the APEX-SZ maps is summarized in 
Table [T] The combination of all contributions 
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to the calibration uncertainty described above re- 
sults in an overall point source flux uncertainty of 
±5.5%. 

4. Data Reduction 

The data consist of 280 bolometer timestreams 
sampled at 100 Hz, telescope pointing data inter- 
polated to the same rate, housekeeping thermom- 
etry data, bolometer bias and readout configu- 
ration data, and other miscellaneous monitoring 
data. The fundamental observation unit is a scan 
comprising twenty 5-s circular sub-scans in AZ /EL 
coordinates, allowing the source to drift through 
the field of view (FOV), as described in §[2] above. 

Data reduction consists of cuts to remove poor- 
quality data, filtering of 1// and correlated noise 
due to atmospheric fluctuations, and binning the 
bolometer data into maps. These steps are de- 
scribed in more detail below. 

4.1. Timestream Data Cuts 

Timestream data are flrst parsed into individual 
circular sub-scans. We reject ^ 7% of the data at 
the beginning and end of the scan where the tele- 
scope deviates from the constant angular veloc- 
ity circular pattern. We reject bolometer channels 
that are optically or electronically unresponsive, 
or lack high-quality flux calibration data; typi- 
cally, 160-200 of the 280 bolometer channels re- 
main after these preliminary cuts. The large num- 
ber of rejected channels is due primarily to low 
fabrication yield for two of the six bolometer sub- 
array wafers. 

We reject spikes and step-like glitches caused 
by cosmic rays or electrical interference. These 
are infrequent and occur on time scales faster than 
the detector optical time constant. We use a sim- 
ple signal-to-noise cut on the data to reject these, 
since the timestream is noise dominated even for 
the ~1 mK Bullet cluster signal. Step- like glitches 
are often correlated across many channels in the 
array, so we reject data from all channels when- 
ever a spike or glitch in > 2% of the channels is 
detected. These cuts result in a loss of 8% of the 
remaining data. 

For each circular sub-scan, we also reject chan- 
nels that have excessive noise in signal band, re- 
sulting in a loss of 19% of the remaining data. 



4.2. Atmospheric Fluctuation Removal 

After the timestream data cuts, fluctuations in 
atmospheric emission produce the dominant sig- 
nal in the raw bolometer timestreams. The atmo- 
spheric signal is highly correlated across the array, 
which can be exploited to remove the signal. Prin- 
cipal component analysis (PCA) has been used by 
some groups to reduce the atmospheric signal (see, 
e.g., Scott et al. 2008; Laurent et al. 2005). How- 
ever, the effect of PCA filtering on the source is 
difficult to predict and is a function of the atmo- 
spheric conditions. We have developed an anal- 
ysis strategy that reduces the atmospheric signal 
through the application of spatial filters that have 
a constant and well understood effect on the sig- 
nals we are attempting to measure. 

Atmospheric fiuctuation power is expected to 
follow a Kolmogorov spatial power spectrum, with 
most power present on scales larger than the sep- 
aration between beams as they pass through the 
atmosphere, resulting in an atmospheric signal 
that is highly correlated across the array. To re- 
duce these fiuctuations, we first remove a poly- 
nomial and an elevation dependent airmass opac- 
ity model from each channel's timestream, then 
remove a first-order two-dimensional spatial poly- 
nomial across the array for each time-step. This 
algorithm is described in detail in the two follow- 
ing subsections. This atmospheric fiuctuation re- 
moval strategy requires that both the spatial ex- 
tent of the scan pattern and the instantaneous ar- 
ray FOV are larger than the source. The 6' radius 
circular scan and the 23' array FOV allow us to 
recover most of the Bullet cluster's fiux, but some 
extended emission is lost as is described in §[5j 

4-.2.1. Timestream Atmosphere Removal 

We observe scan-synchronous signals in the 
bolometer timestreams due to elevation-dependent 
atmospheric emission. The optical path length L 
through the atmosphere is proportional to the 
cosecant of the elevation angle e, L oc csc(e). The 
change in optical path-length is nearly a linear 
function of elevation angle over the 6' circular 
scan radius. In the circular drift scans, this mod- 
ulation of the elevation-dependent opacity pro- 
duces an approximately sinusoidal modulation 
in the bolometer timestream. For each channel 
in each scan, we simultaneously fit and remove 
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Table 1 

APEX-SZ Flux Calibration Uncertainty 



Source 


Uncertainty 


WMAP Mars temperature at 93 GHz 


±1.0% 


Rudy model to WMAP scaling factor at 93 GHz 


±1.0% 


93 GHz to 150 GHz frequency scaling factor 


±0.9% 


Frequency band center 


±1.4% 


Beam solid angle 


±4.0% 


Atmospheric attenuation 


±0.4% 


Gain variation 


±3.0% 


Total 


±5.5% 



an atmospheric model consisting of this cosecant 
function plus an order 20 polynomial (one degree 
of freedom per circular sub-scan) to remove slow 
drifts in the atmospheric opacity and readout 1// 
noise. This scheme effectively removes the com- 
mon scan-induced atmospheric signal as well as 
most of the atmospheric fluctuation power below 
the frequency of the circular sub-scan (0.20 Hz), 
while only modestly affecting the central Bullet 
cluster signal. 

4-2.2. Spatially Correlated Atmosphere Removal 

Removing the cosecant-plus-polynomial model 
from the timestream data reduces low-frequency 
atmospheric fluctuation power, but not higher fre- 
quency power corresponding to smaller spatial 
scales near those where the cluster signal occurs. 
To reduce these fluctuations, the atmosphere can 
be modeled as a spatially correlated signal across 
the array pixel positions on the sky with a low- 
order two-dimensional polynomial function. At 
each time step, we flt and remove a low-order 
two-dimensional polynomial function across the 
axray, similar to the procedure described in Say- 
ers (2007). The relative gain coefficients for each 
bolometer channel are calculated by taking the ra- 
tio of each channel's timestream, which is domi- 
nated by correlated atmospheric noise, to a me- 
dian timestream signal generated from all chan- 
nels. With the favorable atmospheric conditions 
of these observations, we find that a flrst order 
spatial polynomial (offset and tilt) is adequate to 
remove most of the atmospheric signal while pre- 
serving the cluster signal. 



Bolometer channels with excess uncorrelated 
noise are more easily identified after removing the 
correlated atmospheric noise component; we reject 
these noisy channels, then perform the spatially 
correlated signal removal a second time. 

4.3. Map-Making 

The atmospheric removal algorithms described 
above act as a high-pass filter. They suppress sig- 
nals on scales comparable to the scan length or 
the focal plane FOV. The cluster emission can be 
quite extended, and therefore the data reduction 
filtering process attenuates diffuse flux in the clus- 
ter signal and produces small positive sidelobes 
around the cluster decrement. The data reduction 
pipeline filters can be tailored, within limits, to 
meet various scientific objectives. Thus, our pri- 
mary data products consist of two different high 
signal-to-noise maps of the cluster. 

For one map, we mask a circular region cen- 
tered on the cluster source prior to fitting the 
timestream and spatial filters described in § |4.2| 
then apply the resulting filter functions to the en- 
tire data set, including the source region. The 
source-mask procedure prevents the cluster sig- 
nal within the masked region from influencing the 
baseline fits, and thus reduces attenuation of the 
source central decrement and extended emission 
at the expense of increasing the contribution of 
low-frequency noise in the map center. We choose 
a source-mask radius of 4.75' as a compromise 
between attenuation of diffuse emission and in- 
creased map noise. We use the source-masked map 
to visually interpret the morphology and extended 
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emission in the cluster. These results are discussed 
in §|5l] 

We also produce a map in which we do not mask 
the source when applying filters. The non-source- 
masked map is used for model parameter estima- 
tion because it has higher signal to noise in the 
central region of the map. In addition, it is easier 
to take into account the effects of the data reduc- 
tion filters, or transfer function, on the underlying 
sky intensity distribution, which is necessary for 
comparing the data to the model for parameter 
estimation. The fitting procedure and results are 
described in more detail in § |5.2[ 

For each of the two maps, the post-cut, filtered 
timestream data are binned in angular sky coordi- 
nates to create maps. For a given scan, a map 
is created from each bolometer channel, apply- 
ing the channel's pointing offset and flux calibra- 
tion. A coadded scan map is created by combining 
individual channel maps with minimum variance 
weighting in each pixel, using the sample variance 
of the conditioned timestream data in the scan. 
The final coadded map is created by combining all 
scan maps, again with minimum variance weight- 
ing in each pixel. We bin maps at a resolution of 
10" to oversample the beam. The source-masked 
map that we present in §|5]is convolved with a 1' 
FWHM Gaussian to smooth noise fluctuations to 
the angular size of the beam. However, the ra- 
dial profiles presented below and the non-source- 
masked map used for model fitting do not include 
this additional smoothing. 

5. Results 

5.1. Temperature Map 

Figures [2] and [3] show the source-masked tem- 
perature map from our observations of the Bul- 
let cluster. The map has a resolution of 85" 
FWHM which results from the combination of the 
58" instrumental beam, the data reduction filters, 
and a final 1' FWHM Gaussian smoothing of the 
map. The source-masked map is shown in order 
to provide a more accurate representation of the 
extended emission and cluster morphology. The 
noise in the central region of the source-masked 
map is 55 ^K^ms per 85" FWHM resolution ele- 
ment. Near the cluster center, the emission hints 
at elongation in the East- West direction, which is 
along the axis between the main and sub-cluster 



gas detected in the X-ray, see §5.4 The more 



extended emission appears to be elongated in the 
Northwest-Southeast direction, which is the major 
axis of the best-fit elliptical /3 model discussed in 
5.2 Figure [3] shows the centroid position of the 



best-fit elliptical (3 model, and the position of the 
dust obscured, lensed galaxy detected at 270 GHz 
(2008). 



5.3 



by Wilson et al. (2008). As discussed in 
see no evidence for emission from this source in 
our 150 GHz map. 

Radial profiles for the unsmoothed source- 
masked and non-source-masked maps are shown 
in Figure |4] The source-masked map has a signal- 
to-noise of 10 within the central 1' radius, com- 
pared to 23 for the non-source-masked map, due 
to the fact that source-masking allows more large- 
scale atmospheric fluctuation noise to remain in 
the map. However, the source-masked map pre- 
serves signal on larger spatial scales than the non- 
source-masked map. In both the source-masked 
and non- source-masked maps, the sky intensity 
distribution is flltered by the instrument beam 
and data reduction pipeline described in §|4j We 
do not renormalize the map amplitudes, since the 
source is extended and an assumption would need 
to be made about the shape of the sky-brightness 
distribution to do so. However, in order to ac- 
curately estimate cluster parameters such as the 
central temperature decrement, a model for clus- 
ter emission must be adopted, and the instrument 
beam and data reduction filtering must be taken 
into account. 

5.2. Fit to Elliptical j3 Model 

We fit an elliptical (3 model to the non-source- 
masked temperature map to allow a straightfor- 
ward comparison of cluster gas properties derived 
from our measurements to those derived from X- 
ray observations. In all analyses here we assume 
the cluster gas is well-described by an isothermal 
/? model, and is in hydrostatic equilibrium. These 
assumptions are unphysical in the case of the Bul- 
let cluster, which is a dynamically complex merg- 
ing system where the gas is separated from the 
rest of the mass (Clowe et al. 2006) . However, we 
find that with the sensitivity and spatial resolu- 
tion of the observations, these assumptions yield 
an adequate description of the observed emission. 

We model the three-dimensional radial profile 
of the electron density with an isothermal /3 model 



7 



-55 40 
45 



o 50 

o 

o 

CM 

_o 

n3 

.£ 00 
o 

(D 
Q 



05 
10' 




07^^01 



00" 



59" 



58" 



57" 



Right Ascension (J2000) 



-55 40 

45' 



o 50 

o 

o 

CM 

_o 

.£ 00 
o 

(D 
Q 



05 



10 




07'^Or 



00"' 59"' 58"' 57" 
Right Ascension (J2000) 



200 
100 


-100 

-200 I 

-300 =i- 
1-400 
1-500 

-600 



200 
100 


-100 

oa 

-200 I 
-300 ^ 
1-400 

-500 
-600 



Fig. 2. — (Top) Temperature map of the Bullet cluster system from the source-masked data reduction, with 
scale in CMB temperature units. The circle in the lower left corner represents the 85" FWHM map resolution 
which is the result of the instrument beam and data reduction filter convolved with the 1' FWHM Gaussian 
smoothing applied to the map. (Bottom) Difference map made by multiplying alternate scan maps by +1 
and —1, respectively, then coadding all scan maps with minimum variance weighting, in the same manner 
as was used to produce the temperature map shown in the top panel. The contour interval is 100 ^Kcmb in 
both maps. 
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Fig. 3. — Temperature map detail from Figure |2] 
with color scale adjusted to the limits of the detail 
region, and a contour interval of 100 /j,Kcmb- The 
-|- marker indicates the centroid position of the 
best-fit elliptical f3 model, see § |5.2| The * marker 
indicates the position of the bright, dust obscured, 
lensed galaxy detected at 270 GHz by Wilson et al. 
(2008), see 



5.3 



(Cavaliere & Fusco-Femiano 1978): 



2 

ne{r) = rieo ( 1 + ^ 



-3/3/2 



(1) 



Here, rieo is the central electron number density, 
Tc is the core radius of the gas distribution, and /3 
describes the power-law index at large radii. 

The radial surface temperature profile of the 
SZE takes a simple analytic form: 

/ /,2\ (l-3/5)/2 

ATsz = ATo + ^ j ' (2) 

where ATq is the central temperature decrement, 
and 9c — Tc/Da is the core radius divided by the 
angular-diameter distance. A similar form exists 
for the X-ray surface brightness. 

Because of the significant ellipticity in the mea- 
sured SZE intensity profile, we generalize the clus- 
ter gas model to be a spheroidal rather than spher- 
ical function of the spatial coordinates: 



ATsz = ATo(l + A + B) 



(l-3/3)/2 



(3) 



with 



A 



(cos(^)(X - Xq) + sin(^)(y - Yo)Y 



sin(^)(X - Xq) + cos(^)(y - Y^)f 

Mc)2 



Here {X — Xq) and (Y — Yq) are angular offsets 
on the sky in the RA and DEC directions, with 
respect to center positions Xq and Yq. The axial 
ratio, 77, is the ratio between the minor and major 
axis core radii, $ is the angle between the major 
axis and the RA (X) direction. ATq is given by 
the gas pressure integrated along the central line 
of sight through the cluster: 



ATq 

TcMB 



f{x,T,)n,{l)Te{l)dl (4) 



where x — hv/kT, f(x,Te) describes the fre- 
quency dependence of the SZE, ctt is the Thom- 
son scattering cross-section, and Tcmb — 2.728 K. 
For all results in this paper, we use the relativis- 
tic SZE spectrum f{x,Tg) provided by Nozawa 
et al. (2000), and neglect the kinematic effect. At 
150 GHz and Te = 13.9 keV (see §|5l]), this is a 
9% correction to the non-relativistic value. 

To accurately estimate (3 model parameters for 
the cluster, the instrument beam and data reduc- 
tion filters, or transfer function, must be applied 
to the model before comparing it with the data 
(see, e.g., Benson et al. 2003; Reese et al. 2000). 
We characterize this transfer function by creating 
a map from a simulated point source, convolved 
with the instrument beam, and inserted into a 
noiseless timestream, similar to the method de- 
scribed in Scott et al. (2008). The point source 
transfer function map K, is then convolved with a 
simulated /3 model cluster map B to generate a fil- 
tered model map B' , which is a noiseless simulated 
APEX-SZ observation of a /3 model cluster. The 
filtered model map, B' , is then differenced with 
the data map, M , and model parameters are esti- 
mated by minimizing a statistic. 

Simulating maps of many different cluster mod- 
els is required for model parameter fitting. Con- 
volving the cluster model with the point source 
transfer function map is much faster than process- 
ing each model through the reduction pipeline. We 
find that the resulting simulated maps from both 
methods agree sufficiently well to have negligible 
effect on the parameter estimation results. 

We use the unsmoothed non-source-masked 
map with 10" pixelization described in § 4.3 for 



all parameter estimation described below, since 
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this map has lower noise and a more easily char- 
acterized transfer function than the source-masked 
map shown in Figure [2] Diffuse cluster emission is 
more attenuated in the non-source-masked map, 
but this is taken into account using the point 
source transfer function. 

Map noise properties are assessed in the spatial 
frequency domain using jackknife noise maps (see 
Sayers 2007; Sayers et al. 2009). To estimate the 
noise covariance C„, we assume that the noise is 
stationary in the map basis. With this assump- 
tion, the Fourier transform of the noise covariance 
matrix, C„, is diagonal, and the diagonal elements 
are equal to the noise map power spectral den- 
sity (PSD). For each of 500 jackknife noise map 
realizations, we find the two-dimensional Fourier 
transform, then average over all realizations. This 
averaged map PSD is the experimental estimate 
of the diagonal elements of C„. However, these 
jackknife maps do not include fluctuations due to 
the primary CMB anisotropics. We estimate the 
CMB signal covariance from the WMAP5 best- 
fit power spectrum (Nolta et al. 2009) convolved 
with the point source transfer function described 
earlier and add it to the jackknife noise PSD to 
determine the total covariance matrix. 

We construct a statistic for the model fit 
using the transform of the filtered /3 model, B', 
and the transform of the central 14' x 14' portion 
of the data map A4 as: 



(5) 



Using Markov Chain Monte Carlo (MCMC) meth- 
ods, the likelihood, C — e^^x ^ jg sampled in 
the 7-dimensional model parameter space and in- 
tegrated to find the marginal likelihood distribu- 
tions of the P model parameters. The model pa- 
rameter estimates and uncertainties that we report 
are the maximum likelihood values and constant- 
likelihood 68% confidence intervals, respectively, 
of the marginal likelihood distributions. 

The above approach to noise covariance estima- 
tion is chosen for its simplicity and because we do 
not have enough linear combinations of individual 
scan maps to fully sample the noise covariance ma- 
trix using jackknife noise maps. But, the method 
relies on several simplifying assumptions, includ- 
ing that the bolometer noise is stationary for each 
100-s scan, the timestream noise is uncorrelated 
from scan to scan, and the map coverage is uni- 



form. Our map coverage is not actually uniform, 
but we find through simulations of non-uniform 
Gaussian noise maps that the statistic is not 
significantly affected. In addition, the validity of 
the approach is tested by inserting simulated clus- 
ters into the real timestream data; the simulated 
cluster parameters are accurately recovered within 
the estimated uncertainties. 

Results of the (3 model parameter estimation 
are given m Table [2] Due to the de generacy be- 
tween the core radius 9c and /3 parameters, we as- 
sume a prior probabihty density on /3 of 1.04tQ;i[Q, 
which is found from fits to ROSAT X-ray data by 
Ota & Mitsuda (2004). Hallman et al. (2007) find 
that in hydro/N-body simulations, (3 derived from 
fits to SZE profiles is higher than that from X- 
ray, with Psz / Px-my = 1-21 ± 0.13 for fits within 
'^500- We do not account for that factor here due 
to the significant uncertainty in the X-ray derived 
(3 value, but we note that our SZE data prefer a 
higher value for /? than the peak value of the prior. 



We further discuss this choice of prior in § 5.5 



The best-fit (3 model fits the data well, with a 
reduced value of 1.008, and with 7219 degrees of 
freedom (DOF) has a probability to exceed (PTE) 
of 31.5%. The difference map between the data 
map, M, and the best-fit filtered /3 model, B' , 
shows no evidence of residual cluster structure or 
point sources. 

Radial profile plots of the best-fit (3 model, B, 
and the filtered f3 model map, B' , are shown in 
Figure |4] Also plotted for comparison are the 
radially binned data from the unsmoothed non- 
source-masked map A4, used for model fitting, and 
the unsmoothed source-masked map, used to visu- 
alize extended emission (without the 1' Gaussian 
smoothing used in Figures [2] & [3| . Uncertainties 
in both sets of radially binned data are highly cor- 
related due to large-spatial-scale correlated noise 
in the maps. The coincidence of the non-source- 
masked data {M, filled circles) and the filtered 
best-fit (3 model {B' , red solid line) show that the 
data and best-fit /3 model are in good agreement. 
The source-masked map preserves signal on larger 
spatial scales than the non-source-masked map, 
and is useful for visualizing extended emission on 
larger spatial scales. But, as expected, even the 
source-masked map attenuates signal on scales ex- 
ceeding the 4.75' radius of the source masking, and 
thus has a lower signal amplitude when compared 
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Table 2 
P Model Fit Results 



Parameter 


Description 


Value 


Uncertainty" 


^0 


RA centroid position 


06''58'°30.86'' (J2000) 


±7.4" 




DEC centroid position 


-55°56'46.2" (J2000) 


±7.3" 


ATo 


Central temperature decrement 


-771 mKcmb 


±71 A^KcMB 


yo 


Central Comptonization''(re = 13.9 keV) 


3.31 X lO"'^ 


±0.30 X 10-4 


yo 


Central Comptonization''(re = 10.6 keV) 


3.24 X lO^'* 


±0.30 X 10-4 


0c 


Core radius 


142" 


±18" 


V 


Ellipse minor/major core radius ratio 


0.889 


±0.072 


$ 


Ellipse orientation angle 


-52° 


±20° 




Power-law index 


1.15 


±0.13 



'^Quoted uncertainties are 68% confidence intervals in the marginal likelihood distribution for 
each parameter. The uncertainty in ATq includes a statistical uncertainty from the fit of ±57 /iK, 
and a ±5.5% flux calibration uncertainty. The uncertainties in the centroid parameters Xq and Yq 
include a ±4" pointing uncertainty and are given in units of arcseconds on the sky. 

''Central Comptonization, yo, is a derived parameter, assuming an electron temperature of 13.9 
keV from Govoni et al. (2004) and 10.6 keV from Zhang et al. (2006), an SZE observation 
frequency of 152 GHz, and Tcmb = 2.728 K. It is provided to facilitate comparisons with data at 
other wavelengths. 



with the unfiltered /3 model {B, blue dashed line). 

5.3. Radio and IR source contributions 

Radio sources associated with a galaxy clus- 
ter and background IR galaxy sources can have 
a significant impact on the measurement of the 
SZE emission at 150 GHz. We interpret the pub- 
lished results of observations of the Bullet cluster 
at other frequencies and conclude that the mea- 
sured SZE decrement is not significantly contam- 
inated. 

The most important source of potential con- 
fusion is a bright, dust obscured, lensed galaxy 
in the direction of the Bullet cluster recently re- 
ported by Wilson et al. (2008). This source has 
a flux density of 13.5 ± 0.5 mJy at an observ- 
ing frequency of 270 GHz, and is centered at RA 
06'^58"^37.31", DEC -55°57'1.5" (J2000), « 56" 
to the east of the measured SZE centroid position, 
see Figure [3] Assuming a spectral index a = 3, 
where S* oc i^" , we expect a flux density of 1 .94 mJy 
at 150 GHz corresponding to a temperature in- 
crement of AT = 38 /iKcMB in the 1.5 arcmin^ 



APEX-SZ beam solid angle. This lensed source is 
expected to be the dominant contribution to posi- 
tive flux in the direction of the cluster and we have 
repeated the (3 model fit taking it into account. We 
first add the source at its measured position with 
the predicted 150 GHz flux to the SZE f3 model 
and repeat the model fit. As expected, includ- 
ing the point source results in a slightly (~ c/3) 
deeper decrement; however, the of the model 
fit slightly increases. We next allow the flux of the 
point source to vary along with the other model 
parameters and find that values of positive flux are 
a poorer fit than no source at all. Therefore, we 
have no evidence for significant emission from this 
source at 150 GHz. For the results in this paper, 
we use cluster model parameters derived from fits 
that do not include this IR source. 

The Bullet cluster is also associated with a 
number of relatively compact radio sources and 
one of the brightest cluster radio halos yet dis- 
covered. However, these sources are predicted to 
produce negligible temperature increments in the 
APEX-SZ beam when extrapolated to 150 GHz. 
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Fig. 4. — Radial profile of Sunyaev-Zel'dovich ef- 
fect in the Bullet cluster compared to the best- 
fit P model. Points with error bars are the SZE 
data binned in 1' radial bins, from the non-source- 
masked map, {M, filled circles) and the source- 
masked map (open circles). The lines show the 
radial profile of the best-fit /3 model, unfiltered 
{B, blue dashed line) and after convolving with the 
instrument beam and non-source-masked data re- 
duction filters {B' , red solid line). The non-source- 
masked map radial profile is reasonably well fit 
by the filtered P model. The source-masked map 
preserves signal on larger spatial scales than the 
non-source-masked map, but still attenuates sig- 
nal on scales exceeding the 4.75' radius of the 
source masking. The source-masked data thus 
have a lower signal amplitude when compared with 
the unfiltered (3 model, as expected. See text for 
details. 

Liang et al. (2000) report the detection of 
eight radio point sources all of which have steeply 
falling spectra. Only two of these sources were 
detected with ACTA at 8.8 GHz, and they were 
found to have fiux densities of 3.2 ± 0.5 mJy and 
3.3 ± 0.5 mJy. The spectra of these sources, mea- 
sured between 4.9 and 8.8 GHz are falling with 
a ~ —0.93 and a — —1.33, respectively. Ex- 
trapolating to 150 GHz, the fiux of these sources 
are expected to be 0.24 and 0.07 mJy, correspond- 
ing to CMB temperature increments of AT — 4.7 
and 1.4 //Kcmb in the 1.5 arcmin^ APEX-SZ beam 
soUd angle. 

The radio halo in the Bullet cluster is very lu- 
minous, but has a characteristically steeply falling 



spectrum. Liang et al. (2000) measure the flux 
and spectra for the two main spatial components 
of the halo. At 8.8 GHz, they find the two com- 
ponents to have fiuxes of 3.5 and 0.55 mJy, with 
spectral indices of a = —1.3 and a = —1.4, re- 
spectively. Extrapolating to 150 GHz, the com- 
bined flux from the radio halo is expected to 
be ~ 0.1 mJy. This emission is spread over an 
area comparable to the size of the cluster and 
therefore corresponds to a temperature increment 
< 1 /iKcMB in the APEX-SZ beam. 

5.4. Comparison with X-ray Data 

X-ray emission in the ionized intracluster gas is 
dominated by thermal bremsstrahlung. The X-ray 
surface brightness can be written 

Sx - ^^(/^^^4 Jnen^Ae^dl, (6) 

where n^^i are the electron and ion densities in 
this gas, Aei is the X-ray cooling function, and 
the integral is taken along the line of sight. The 
X-ray flux is proportional to the line-of-sight inte- 
gral of the square of the electron density, resulting 
in emission that is more sensitive to local density 
concentrations than the SZE emission. 

The Bullet sub-cluster and the bow shock are 
apparent in the X-ray surface brightness map 
shown in Figure |5] The SZE contour map of the 
Bullet cluster in Figures |2] & |3] is overlaid on an 
X-ray map and weak lensing surface mass density 
reconstruction from Clowe et al. (2006) |^ The 
X-ray map is made from XMM data (observation 
Id: 0112980201) extracted in the [0.5-2] keV band, 
corresponding to Bullet rest frame energies where 
the X-ray cooling function for hot gas is relatively 
insensitive to temperature. The map is smoothed 
with a 12" Gaussian kernel. 

The SZE contours do not resolve the sub- 
cluster. However, an elongation of the inner con- 
tours to the West suggests that a contribution 
from it may be detected. The observed SZE map 
is consistent with expectations, given the 85" res- 
olution of the SZE map, the different dependence 
of the X-ray and SZE signals on gas density, and 
the mass and temperature difference between the 
two merging components which predict a factor 

^Data are publicly available at 

|http://flamingos.astro.ufl.edu/le0657/public.html[ 
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of ^ 10 lower integrated pressure from the sub- 
cluster. There is no evidence in the SZE contours 
of a contribution from the lensed sub-mm bright 
galaxy discussed in §|5.3| 



-55 52 



HRI (~ 2" resolution) measurements of the in- 
ner 6' radius of the Bullet cluster: (3 




Right Ascension (J2000) 

Fig. 5. — The SZE map of the Bullet system from 
this work, in white contours, overlaid on an X- 
ray map from XMM observations. The green con- 
tours show the weak lensing surface mass density 
reconstruction from Clowe et al. (2006) . The SZE 
contour interval is 100 /iKcMB- 



5.5. Mass- Weighted Temperature 

The combination of cluster SZE and X-ray mea- 
surements can be used to place constraints on the 
thermal structure of the intracluster gas. The SZE 
intensity is proportional to the product of the elec- 
tron density and the electron temperature along 
the line of sight, see equation Q. Therefore, if 
the electron density is known from another mea- 
surement, the SZE can be used to measure a mass- 
weighted temperature. For simplicity, we assume 
here that the intracluster gas is isothermal, but 
a more detailed comparison of the SZE surface 
brightness and projected density could be used to 
constrain the thermal structure in the cluster. 

We perform this calculation with two differ- 
ent descriptions of the intracluster gas density. 
First, we model the spatial distribution of the 
intracluster g spherical /3 model follow- 

ing equation ([ij. We use (3 model parameters 
from Ota & Mitsuda (2004), derived from ROSAT 



— 110 c; + 15.6// 

— iiZ.0_;^Q4 , 



HO. 3 



1 04+''-is 

^•^^-0.10' 



and TieO '^■^^Q 3 ^ -LVJ ^LLL 

We construct an X-ray derived SZE surface 
brightness model from the j3 model electron sur- 
face density profile using equation ^ . To account 
for (3 model uncertainties, we incorporate the val- 
ues and uncertainties for (3, 9c, and rieo as inde- 
pendent Gaussian prio rs. W e then use the analysis 

on the 



method described in S5.2 



to minimize x 



difference between the X-ray derived SZE model, 
convolved with the point source transfer function, 
and the APEX-SZ non-source-masked data. The 
free parameters in the fit are the three (3 model 
parameters, the mass- weighted electron tempera- 
ture Tjng, and the relative map alignment in RA 
and DEC. 

We find Tmg = 11.4 ± 1.4 keV after marginal- 
izing over the other parameters in the fit and 
including the SZE flux calibration uncertainty 
and the effect of the APEX-SZ band center fre- 
quency uncertainty on the relativistic SZE spec- 
trum f{x,Tmg)- The reduced of the best fit 
model is 1.008 with an associated PTE of 31.3%, 
indicating that the spherical (3 density model and 
the assumption of isothermality produce an ac- 
ceptable fit to the data. 

Given the complex morphology of this merging 
system, the validity of the spherical /? model is 
questionable. We therefore repeat the determina- 
tion of the mass-weighted temperature by directly 
comparing X-ray measurements of the projected 
intracluster gas density with the measured SZE 
signal in order to produce a less model-dependent 
measurement of the mass-weighted temperature. 
We make use of the publicly availably electron 
surface density map derived from Chandra X-ray 
satellite data presented in Clowe et al. (2006). 
Using the same analysis as above, and marginal- 
izing over the relative map alignment parame- 
ters, we find mass- weighted electron temperature 
= 10.8 ± 0.9 keV. The fit to the data is again 
good, with a reduced — 1.037 and a PTE of 
20.6%. This is in excellent agreement with the 
value Tmg = 11.4 ± 1.4keV found from the above 
/3 model analysis. 

Given the complex dynamics in the Bullet clus- 
ter, there have been several studies of the tem- 

"^http: / /flamingos. astro. ufl.edu/le0657/public. html 
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perature structure (e.g., Finoguenov, Bohringer, 
& Zhang 2005; Markevitch 2006; Andersson, Pe- 
terson, & Madejski 2007). There have also been 
several pubhshed results for the spectroscopic tem- 
perature within annuli about the cluster center. 
Chandra data was used by Govoni et al. (2004) 
to determine a spectroscopic X-ray temperature 
of Tspec = 13.9 ± 0.7 keV within 0.75 Mpc of 
the cluster center. From the analysis of XMM 
data within an annulus of 0.14 — 0.7 Mpc ra- 
dius, Zhang et al. (2006) find a temperature of 
Tspec = 10.6 ± 0.2 keV. Analyzing the combina- 
tion of XMM and RXTE data, Petrosian, Made- 
jski, & LuH (2006), find T.p^c = 12.1 ± 0.2 keV 
within a radius of 0.95 Mpc. The published X- 
ray spectroscopic temperatures span a range much 
larger than the stated uncertainties in the mea- 
surements. Given the complex thermal structure 
for the cluster, and the presence of gas at tem- 
peratures corresponding to energies at or above 
the upper limits of the Chandra and XMM en- 
ergy response, the variation in the measured X-ray 
spectroscopic temperature is not surprising. The 
mass-weighted temperature found with APEX-SZ 
falls near the lowest of the reported X-ray spec- 
troscopic temperatures. However, we do not ex- 
pect exact agreement between the mass-weighted 
and spectroscopic temperatures. Using Chandra 
data for a sample of 13 relaxed clusters, Vikhlinin 
et al. (2006) find that, due to the presence of 
thermal structure in the intracluster gas, the X- 
ray spectroscopic temperature is typically a fac- 
tor of Tspec/Tmg = l-H ± 0.06 larger than the 
X-ray derived mass-weighted electron tempera- 
ture. This is consistent with the simulation re- 
sults of Nagai, Vikhhnin, & Kravtsov (2007) who 
find T^rier/T.nn ~ 1-14 for relaxed clusters and 



spec 

Tspec/Tmg ~ 1.12, with a somewhat larger scatter, 
for unrelaxed systems. Naively applying this cor- 
rection to the published X-ray spectroscopic tem- 
peratures, we infer results for mass-weighted tem- 
peratures that bracket the APEX-SZ result. 

5.6. Gas Mass Fraction 

Using the SZE measurements, we construct a 
model for the intracluster gas distribution which, 
when combined with X-ray measurements, can be 
used to determine the gas mass, total mass, and 
therefore gas mass fraction of the cluster. The 
gas mass is estimated by integrating a spheroidal 



model for the cluster gas, following LaRoque et al. 
(2006). 

Several assumptions about the model must be 
made to estimate the gas mass. We assume that 
the cluster gas is isothermal in order to convert 
pressure to density. We also assume spheroidal 
symmetry for the gas distribution in order to con- 
vert the two-dimensional SZE integrated pressure 
measurement to a three-dimensional gas distribu- 
tion. 

We consider two simple cases, an oblate spheroid 
generated by rotation about the minor axis and a 
prolate spheroid generated by rotation about the 
major axis, where the symmetry axis is in the sky 
plane. The gas mass, under these assumptions, 
becomes: 



^gas(f) = S^ieneO'mpDA^ 



t/Da 



dX dY dZ 

2\ -3/3/2 




(7) 



where /ie is the nuclcon/clcctron ratio, taken to 
be 1.16 (Grego et al. 2001). The factor of eight 
is due to integrating over only one octant of the 
spheroid. The factor ( is set to unity in the case 
of oblate spheroidal symmetry, while in the case 
of prolate spheroidal symmetry, C is set to ry. 

The total cluster mass is estimated by assum- 
ing hydrostatic equilibrium and integrating the in- 
ferred gas distribution (Grego et al. 2000) to find: 



Ptotal — 



fcTe 

47rG/imp 



V^lnpg 



(8) 



Here /x is the mean molecular weight of the intr- 
acluster gas, which is assumed to be 0.62 (Zhang 
et al. 2006). Using equations ([t]) & (|8|, and our 
model parameters in Table [2] we calculate the gas 
mass, total mass, and gas mass fraction for the 
cluster. In Table [3] we give these results. The gas 
mass fraction results for a prolate gas distribution 
model are ~ 3% larger than those for an oblate 
model, while the total mass and gas mass are 
< 18% larger. We quote only the oblate spheroidal 
results. 

We calculate our results within two different 
radii. The first is the radius of the cluster at 
which its mean density is equal to 2500 times 
the critical density at the redshift of the cluster. 
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''2500 ■ The second radius is 1.42 Mpc, which is 
the same radius used by Zhang et al. (2006) for 
their gas mass fraction calculation. This will al- 
low for a more direct comparison to their result, 
and is also near where our measured SZE radial 
profile has unity signal to noise. For all results, 
we assume a ACDM cosmology, with h = 0.7, 
^lyn = 0.27, and f^A = 0.73. The results of this 
analysis are summarized in Table [3] Under the 
assumption of a 10.6 keV mass- weighted temper- 
ature (the lowest of the published X-ray spectro- 
scopic temperatures and near our mass-weighted 
temperature results in 



fractions fg 



5.51 



we find gas mass 
and 0.179 ± 0.036 



0.216 ± 0.031 
within r25oo and 1.42 Mpc, respectively. The fact 
that the computed gas fraction in the central re- 
gion significantly exceeds the cosmic average de- 
termined by WMAP5 {fg = 0.165 ± 0.009, Dunk- 
ley et al. 2009), and a lower value observed in 
relaxed clusters {fg ~ 0.12, see, e.g., VikhUnin 
et al. 2009), is Ukely due to deviations of the in- 
tracluster gas from isothermal hydrostatic equi- 
librium. On larger scales, baryon fractions pro- 
duced for the range of reported X-ray tempera- 
tures bracket the pubhshed results using X-ray 
and weak lensing data. Bradac et al. (2006) 
measure a gas mass fraction fg = 0.14 ± 0.03 by 
comparing the gas mass calculated from Chandra 
X-ray measurements to weak lensing total mass 
measurements in a 4.9' x 3.2' box roughly cen- 
tered around the cluster. Zhang et al. (2006), 
measured a gas mass fraction fg = 0.161 ± 0.018 
within a radius of 1.42 Mpc. Despite the limita- 
tions of applying a hydrostatic equilibrium model 
to this merging cluster, the APEX-SZ results for 
the gas mass fraction are in good agreement with 
previous work. 

6. Conclusions 

Measurements of the SZE provide a robust and 
independent probe of the intracluster gas proper- 
ties in galaxy clusters. The APEX-SZ 150 GHz 
observations detect the Bullet system with 23 cr 
significance within the central 1' radius of the 
SZE centroid position. We do not expect to see 
a resolved signal from the Bullet sub-cluster in 
the 150 GHz 85" FWHM resolution SZE maps, 
and no obvious feature, such as a secondary peak, 
is present. We expect no significant contamina- 
tion of the observed SZE decrement due to radio 



sources, and there is no evidence for significant 
contamination by a known bright lensed dusty 
galaxy. 

We process an elliptical /? model through the 
observation transfer function and fit it to the mea- 
sured temperature decrement map. We also mea- 
sure the cluster mass-weighted electron temper- 
ature and gas mass fraction with the SZE data. 
Combining the APEX-SZ map with a map of pro- 
jected electron surface density from Chandra X- 
ray observations, we determine the mass- weighted 
temperature of the cluster gas to be Tmg — 10.8 ± 
0.9 keV. This value is consistent with the lowest 
X-ray spectroscopic temperatures reported for this 
cluster and should be less sensitive to the details of 
the cluster thermal structure. The derived baryon 
fraction is also found to be in reasonable agree- 
ment with previous X-ray and weak lensing deter- 
minations. 

Throughout this work, we make an assump- 
tion of isothermal cluster gas. Clearly, incorpo- 
rating thermal structure, measured by X-ray ob- 
servations, in the analysis of the SZE data would 
improve the determination of the gas distribution 
and gas mass fraction. Ultimately, a more sophis- 
ticated analysis could be implemented that com- 
bines X-ray, SZE, and weak lensing data and re- 
laxes assumptions of hydrostatic equilibrium be- 
tween the gas and dark matter components of the 
cluster. This is particularly important for a de- 
tailed understanding of actively merging systems 
such as the Bullet cluster. 
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Table 3 

Mass Estimates for the Bullet Cluster System 



Te(keV)^ Mean Overdensity T-jnt (') T-int (l^Pc)^ Gas Mass Fraction Gas Mass Total Mass (IO^^Mq) 

13.9 ± 0.7 2506 ± 233 2.77 0.739 0.124 ± 0.022 0.944 ± 0.105 7.56 ± 0.70 

13.9 ± 0.7 961 ± 98 5.32 1.42 0.106 ± 0.024 2.20 ± 0.33 20.6 ±2.1 

10.6 ± 0.2 2521 ± 230 2.15 0.572 0.216 ± 0.031 0.765 ± 0.072 3.54 ± 0.32 

10.6 ± 0.2 734 ± 66 5.32 1.42 0.179 ± 0.036 2.83 ± 0.40 15.7 ±1.4 



Note. — Two different isothermal electron temperatures are assumed, in order to bracket the range of X-ray spectroscopic temperatures reported 
in the literature. The top two rows assume an isothermal electron temperature of 13.9 ± 0.7 keV. The bottom two rows assume an isothermal 
electron temperature of 10.6 ± 0.2 keV. For each electron temperature, we integrate to ^^2500' radius within which the mean cluster density 
is 2500 times greater than the critical density at the redshift of the cluster. For each electron temperature, we also integrate to a fixed radius 
of 1.42 Mpc, allowing a direct comparison to results in Zhang et al. (2006). This radius is also near where our measured SZE radial profile has 
unity signal to noise. For all results, we assume a ACDM cosmology, with h = 0.7, f^m = 0.27, and Ctj^ = 0.73. Uncertainties in the gas mass 
and gas mass fraction include a ±5.5% SZE flux calibration uncertainty. 

^Isothermal electron temperature 

'^Angular integration radius. 

''Physical integration radius. 
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